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Abstract 



The contrast between model predictions for the transverse momentum spectra of 
J/ip observed in Au-Au collisions at RHIC is extended to include effects of nuclear 
absorption. We find that the difference between initial production and recombina- 
tion is enhanced in the most central collisions. Models utilizing a combination of 
these sources may eventually be able to place constraints on their relative magni- 
tudes. 
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The original motivation [1] to consider heavy quarkonium production in a 
region of color deconfinement was the realization that the number of charm 
pairs initially produced in a central heavy collision will be at least of order 10 
[2,3] at RHIC 200 energy. Preliminary model estimates indicated that interac- 
tions of a quark and antiquark which were not produced in the same nucleon- 
nucleon collision (the so-called off-diagonal pairs) will result in formation of 
additional J/ip (now called recombination) which would mask the expected 
suppression [4] due to medium effects. However, these model estimates proved 
to be very sensitive to details of various parameters which control the recom- 
bination process [5]. PHENIX results have now provided a first look at the 
total J I ip population in terms of the ratio Raa over a large range of collision 
centrality [6] . Surprisingly, one sees a suppression factor which is very close to 
that observed at SPS energy. A number of models predicted that there would 
be substantially more suppression at RHIC, and all were thus incompatible 
with these initial results. However, models in which some sort of in-medium 
recombination process is included find a total J/ip population which is not 
incompatible with the PHENIX results. This can be taken as evidence that 
the recombination process is indeed significant at RHIC energy. However, one 
must have confidence that the suppression parameters used are indeed ap- 
propriate at RHIC, which at this time has no independent confirmation. In 
fact, there is a proposed alternate scenario [7,8] in which the observed sup- 
pression at both SPS and RHIC energies is entirely due to dissociation of the 
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ip' and Xc which would otherwise have provided additional J/i/j from decay 
feed-down. We show in Figure 1 some model predictions with and without 
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Fig. 1. Comparison of dissociation and recombination model predictions for 
J/tp R A a ■ 



recombination compared with the initial PHENIX measurements. The set of 
model parameters was motivated in [5] which uses an initial temperature T Q 
= 500 MeV and N c5 = 10 initially-produced cc pairs for central collisions, 
where the charm quark momentum distribution is calculated in pQCD. The 
prediction for no recombination (stars) leads to almost complete suppression 
of initially-produced J/ip . Inclusion of in-medium charm quark recombina- 
tion (blue circles) dominates over almost the complete centrality range, and 
the sum (green diamonds) is seen to follow the trend of the measured points 
fairly well. We note that since the normalization of Raa was not contained 
in the model predictions, we choose to normalize Raa to unity at the most 
peripheral points. 

To illustrate the fragility of such an interpretation, we also show model pre- 
dictions in which no recombination process occurs and the initial temperature 
is decreased to 350 MeV. The resulting Raa (black squares) is also seen to be 
roughly compatible with the data. Thus in the absence of independent infor- 
mation on the initial temperature parameter in this model, we cannot use the 
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Raa measurements to make a definitive statement. 

Fortunatately, the transverse momentum spectra of the surviving J /if} are 
also measured, and we can use this additional information to test the var- 
ious scenarios. In the next section we review the use of spectra as a signal 
of the production mechanism. The following section summarizes the crucial 
role of centrality dependence and presents a comparison with the PHENIX 
measurements. Finally, we include for the first time the effects of nuclear ab- 
sorption on these signals, and illustrate the experimental precision necessary 
for a definitive result. 



2 Spectra as signals of recombination 



The number of J /if) formed in a medium via recombination of pairs of c and 
c quarks is determined by the competition between the recombination and 
dissociation cross sections. In the case where the total number Nj/^ is much 
less than the initial number of quark pairs N c5 , the solution is of the form 

*/ 

Nj/*(t f ) = e{t s )[N m {t Q ) + NlJ A F [V(t) e^)]" 1 dt], (1) 

to 



where to an d tf define the lifetime of the deconfined region, and Af,d are the 
formation and dissociation cross section reactivities, respectively, and V(t) is 
the (expending) volume of the deconfinement region. Note that the function 

e(t f ) = e ~^o XDPsdt would be the suppression factor in this scenario if the 
formation mechanism were neglected. 

One sees that the second term, quadratic in N c5 , is precisely the total num- 
ber of possible recombinations which could occur in the deconfinement vol- 
ume, modified by the factor e(t/)/e(t), which is just the suppression factor for 
J/ip formed between times t and tf. To calculate the transverse momentum or 
rapidity spectra of the J/ip formed via recombination, one must use a version 
of this solution which is differential in the J /if} momentum. Since the cross 
section used here is simple gluo-ionization of the J/ if) bound state by the inter- 
action of its color dipole with free gluons in the medium, the differential cross 
section expression is trivial. To solve the equation, we first generate a set of 
cc pairs according to pQCD amplitudes [9] to fix the kinematic distributions 
of the recombining quarks. The resulting spectra of the J /if} is then 

dNj^ f dt da 
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where the sum over all cc pairs incorporates the total formation reactivity. 

The charm quark spectra we calculate from pQCD have used collinear par- 
ton interactions only. To simulate the effects of confinement in nucleons, and 
more importantly the initial state broadening due to multiple scattering 
of nucleons, we have supplemented the quark pair pt distributions by adding 
a transverse momentum kick to each quark in a diagonal pair, chosen from a 
Gaussian distribution with width characterized by (/ct 2 )- This effect is mag- 
nified in the difference between diagonal and off-diagonal pairs. Since the az- 
imuthal direction of the kx is uncorrelated from pair-to-pair, an off-diagonal 
pair with initial (pt 2 ) will inherit an increase of 2 (kr 2 )- However, the az- 
imuthal correlation inherent in the diagonal pair will result in an increase 
4 (kx 2 )- This effect is evident when we investigate the (&t 2 ) dependence of 
charm quark pairs and J/ip formed from recombination. We find that: 

(PT 2 ) diagmal ccbar = 0.8 + 4 (k T 2 ) GeV 2 
(PT 2 ) off ^ agonal ccbar = 4.9 + 2 (k T 2 ) GeV 2 

(PT )in~medium J/V formation = 2.4 + (k T ) GeV 

We use reference data on J ftp production in pp interactions [10], which of 
course involves diagonal pairs, to fix (kT 2 ) pp = 0.5 ±0.1 GeV 2 . The corre- 
sponding value in nuclear interactions depends on the centrality, which we 
determine in the next section. 



3 Centrality Dependence of (pr 2 ) 

We use the standard random walk picture of initial state NN interactions to 
model the transverse momentum broadening in pA and AA interactions [11]. 
This uses a parameter A 2 , which is proportional to the average (kx 2 ) per initial 
state collision (a factor of 0.5 times the probability for initial state parton 
elastic scattering in a given inelastic NN collision is conventionally absorbed 
in the definition). The net broadening in a pA collision can be parameterized 
in terms of the product of this A 2 and N n , the average number of inelastic 
collisions each nucleon incurs before the interaction in which the charm quark 
pair is produced. In terms of the total number of collisions per initial state 
nucleon (n), this number is 

"n = ±X>-l) = ^. (3) 
n m=l Z 

We extract A 2 = 0.56 + 0.08 GeV 2 from the nuclear broadening measurements 
of J/ip produced in minimum bias d-Au interactions [10]. The corresponding 
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{^t 2 )aa is then calculated as a function of centrality, parameterized by the 
total number N co u of binary collisions. The final results for (pt 2 ) of initially- 
produced J/ip from diagonal cc pairs and (pt 2 ) of J/ip formed in-medium 
are shown in Figure 2. One sees that they differ by increasing amounts as 
centrality increases, which is a direct consequence of their different cc pair 
precursors (diagonal or off-diagonal). We note that the initial production curve 
(higher (pt 2 )) should describe the sequential production scenario, in which all 
observed J/ ip originate from unsupressed direct production [8] . The in- medium 
production curve (lower (pt 2 )) would describe a situation in which all of the 
initially-produced J/ ip have been screened away, and thus all observed J / tp are 
due to in-medium formation processes. We should stress that these calculations 
of (pr 2 ) use the charm quark momentum distributions which follow from initial 
production via pQCD amplitudes. One might expect that the actual physical 
situation would be some weighted sum of these two extreme cases. (For a 
discussion of comparison with preliminary data, see [12].) We should stress 
that these calculations of (pr 2 ) use the charm quark momentum distributions 
which follow from initial production via pQCD amplitudes. 




Fig. 2. Centrality dependence of (pt 2 ) comparing model predictions for initial (di- 
rect) production and in-medium formation, both assuming charm quarks with mo- 
mentum spectra as predicted by pQCD. 
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4 Effects of Nuclear Absorption 

There is an additional effect on the p T spectrum of directly-produced J / ip due 
to their interaction in cold nuclear matter which they encounter on their path 
out of the interaction region. We note that this effect cannot modify the 
J/ip formed via recombination in the medium, since the formation process 
takes place at a later stage. One could include the nuclear absorption effect in 
a full Glauber calculation of the initial production of J/ip , but for simplicity 
we here extract a correction factor for the pt spectra using an effective absorp- 
tion cross section in a path length geometry. This will involve an extension of 
the random-walk picture of individual NN collisions which occur both before 
and after the point at which a hard interaction produces the cc pair which 
eventually emerges as a J ftp. What we wish to calculate in this picture is the 
average number of such NN collisions which occur before the collision in which 
the J j ip production is initiated. In the absence of nuclear absorption, the rela- 
tive probability of observing a J/ip will be independent of the spatial position 
of its production point. This leads to the expression in Eq. 3. With nuclear 
absorption, however, one must weight each of the possible production points 
by the probability that the J/ip produced at a given position will survive to 
be observed in the final state. 

Let P mn be the relative probability that the J/ip survives its path through 
the remaining nuclear matter, where n is the total number of NN interactions 
and m specifies the number in this sequence at which the J/ip production 
point occurs. We parameterize this in a simple straight-line trajectory model, 
to obtain 

Pmn = £ n n ~ m , where e n = exp(—paL rnax /n) is the average absorption factor 
encountered by the J/ip in each path length between NN collisions. Here p is 
the density of nuclear matter and a is the absorption cross section. 

Then the average effective number of NN collisions which an observed J/ip 
will have experienced is 

1 n 

Nn = -^Y,( m - l ) P ™m (4) 

r to t m= i 

where 

Ptot = X/ p>mn = 7~- (5) 

m=l 6 n ~ 1 

In the limiting case of zero nuclear absorption, P mn = 1, P tot = n, and one 
recovers N n = In the opposite limit of infinite absorption, P mn vanishes 
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unless m=n, resulting in N n = n-1. This correctly describes the effective num- 
ber of initial state collisions for a J /if) which is produced in the final collision 
point at the back of the nucleus, the only possibility which avoids complete 
absorption. We show in Figure 3 the corresponding (maximum) increase in 

(PT 2 ). 




Fig. 3. Maximum increase of {pt 2 } f° r initial (direct) production in the limit of 
infinite nuclear absorption. 

For finite nuclear absorption in terms of /? = p a L max , a numerical evaluation 
yields N n » (1 + 0.2 (3) *=±, for < (3 < 1.5 . We show in Figure 4 the 
predicted (pr 2 ) spectra of initially produced diagonal cc pairs, as appropriate 
for initial production of J /if) in Au-Au interactions. One sees that the effect 
of nuclear absorption always results in increased (pt 2 ), which enhances the 
difference with respect to in-medium J /if) production. Even the smallest ab- 
sorption cross section considered (3 mb) shifts the prediction upward to an 
extent which may prove difficult to reconcile with preliminary PHENIX data 
at both forward and mid rapidity [12]. 
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Fig. 4. Variation with nuclear absorption of (pt 2 ) for initial (direct) production. 
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